Search for: All records

Abstract Observations of GeV gamma-ray emission from the well-studied mixed-morphology supernova remnant (SNR) W44 by Fermi-Large Area Telescope and AGILE imply that it is a site of significant cosmic-ray acceleration. The spectral energy distribution (SED) derived from the GeV data suggests that the gamma-ray emission likely originates from the decay of neutral pions generated by cosmic-ray interactions. It is essential to measure the SED of W44 in the X-ray and very-high-energy (VHE) gamma-ray bands to verify the hadronic origin of the emission and to gauge the potential contributions from leptonic emission. We report an upper limit of the nonthermal X-ray flux from W44 of 5  × 10⁻¹³erg cm⁻²s⁻¹in the 0.5–8.0 keV band based on  ∼300 ks of XMM-Newton observations. The X-ray upper limit is consistent with previously estimated hadronic models, but in tension with the leptonic models. We estimate the VHE flux upper limit of  ∼1.2  × 10⁻¹²erg s⁻¹cm⁻²in the 0.5–5.0 TeV range from W44 using data from the Very Energetic Radiation Imaging Telescope Array System. Our nondetection of W44 at VHE wavelengths is in agreement with observations from other imaging atmospheric Cherenkov telescopes and is perhaps consistent with the evolutionary stage of the SNR. 

Abstract Pulsar halos are regions around middle-aged pulsars extending out to tens of parsecs. The large extent of the halos and well-defined central cosmic-ray accelerators make this new class of Galactic sources an ideal laboratory for studying cosmic-ray transport. LHAASO J0621+3755 is a candidate pulsar halo associated with the middle-aged gamma-ray pulsar PSR J0622+3749. We observed LHAASO J0621+3755 with VERITAS and XMM-Newton in the TeV and X-ray bands, respectively. For this work, we developed a novel background estimation technique for imaging atmospheric Cherenkov telescope observations of such extended sources. No halo emission was detected with VERITAS (0.3–10 TeV) or XMM-Newton (2–7 keV) within 1^∘and $10^{\prime }$around PSR J0622+3749, respectively. Combined with the LHAASO Kilometer Square Array (KM2A) and Fermi-LAT data, VERITAS flux upper limits establish a spectral break at  ∼1–10 TeV, a unique feature compared with Geminga, the most studied pulsar halo. We model the gamma-ray spectrum and LHAASO-KM2A surface brightness as inverse Compton emission and find suppressed diffusion around the pulsar, similar to Geminga. A smaller diffusion suppression zone and harder electron injection spectrum than Geminga are necessary to reproduce the spectral cutoff. A magnetic field ≤1μG is required by our XMM-Newton observation and synchrotron spectral modeling, consistent with Geminga. Our findings support slower diffusion and lower magnetic field around pulsar halos than the Galactic averages, hinting at magnetohydrodynamic turbulence around pulsars. Additionally, we report the detection of an X-ray point source spatially coincident with PSR J0622+3749, whose periodicity is consistent with the gamma-ray spin period of 333.2 ms. The soft spectrum of this source suggests a thermal origin. 

Abstract G106.3+2.7, commonly considered to be a composite supernova remnant (SNR), is characterized by a boomerang-shaped pulsar wind nebula (PWN) and two distinct (“head” and “tail”) regions in the radio band. A discovery of very-high-energy gamma-ray emission (E_γ> 100 GeV) followed by the recent detection of ultrahigh-energy gamma-ray emission (E_γ> 100 TeV) from the tail region suggests that G106.3+2.7 is a PeVatron candidate. We present a comprehensive multiwavelength study of the Boomerang PWN (100″ around PSR J2229+6114) using archival radio and Chandra data obtained two decades ago, a new NuSTAR X-ray observation from 2020, and upper limits on gamma-ray fluxes obtained by Fermi-LAT and VERITAS observatories. The NuSTAR observation allowed us to detect a 51.67 ms spin period from the pulsar PSR J2229+6114 and the PWN emission characterized by a power-law model with Γ = 1.52 ± 0.06 up to 20 keV. Contrary to the previous radio study by Kothes et al., we prefer a much lower PWNB-field (B∼ 3μG) and larger distance (d∼ 8 kpc) based on (1) the nonvarying X-ray flux over the last two decades, (2) the energy-dependent X-ray size of the PWN resulting from synchrotron burn-off, and (3) the multiwavelength spectral energy distribution (SED) data. Our SED model suggests that the PWN is currently re-expanding after being compressed by the SNR reverse shock ∼1000 yr ago. In this case, the head region should be formed by GeV–TeV electrons injected earlier by the pulsar propagating into the low-density environment. 

Abstract Dark matter is a key piece of the current cosmological scenario, with weakly interacting massive particles (WIMPs) a leading dark matter candidate. WIMPs have not been detected in their conventional parameter space (100 GeV ≲ M χ ≲ 100 TeV), a mass range accessible with current Imaging Atmospheric Cherenkov Telescopes. As ultraheavy dark matter (UHDM; M χ ≳ 100 TeV) has been suggested as an underexplored alternative to the WIMP paradigm, we search for an indirect dark matter annihilation signal in a higher mass range (up to 30 PeV) with the VERITAS γ -ray observatory. With 216 hr of observations of four dwarf spheroidal galaxies, we perform an unbinned likelihood analysis. We find no evidence of a γ -ray signal from UHDM annihilation above the background fluctuation for any individual dwarf galaxy nor for a joint-fit analysis, and consequently constrain the velocity-weighted annihilation cross section of UHDM for dark matter particle masses between 1 TeV and 30 PeV. We additionally set constraints on the allowed radius of a composite UHDM particle. 

Context. The response of imaging atmospheric Cherenkov telescopes to incident γ -ray-initiated showers in the atmosphere changes as the telescopes age due to exposure to light and weather. These aging processes affect the reconstructed energies of the events and γ -ray fluxes. Aims. This work discusses the implementation of signal calibration methods for the Very Energetic Radiation Imaging Telescope Array System (VERITAS) to account for changes in the optical throughput and detector performance over time. Methods. The total throughput of a Cherenkov telescope is the product of camera-dependent factors, such as the photomultiplier tube gains and their quantum efficiencies, and the mirror reflectivity and Winston cone response to incoming radiation. This document summarizes different methods to determine how the camera gains and mirror reflectivity have evolved over time and how we can calibrate this changing throughput in reconstruction pipelines for imaging atmospheric Cherenkov telescopes. The implementation is validated against seven years of observations with the VERITAS telescopes of the Crab Nebula, which is a reference object in very-high-energy astronomy. Results. Regular optical throughput monitoring and the corresponding signal calibrations are found to be critical for the reconstruction of extensive air shower images. The proposed implementation is applied as a correction to the signals of the photomultiplier tubes in the telescope simulation to produce fine-tuned instrument response functions. This method is shown to be effective for calibrating the acquired γ -ray data and for recovering the correct energy of the events and photon fluxes. At the same time, it keeps the computational effort of generating Monte Carlo simulations for instrument response functions affordably low. 

Abstract The Breakthrough Listen Initiative is conducting a program using multiple telescopes around the world to search for “technosignatures”: artificial transmitters of extraterrestrial origin from beyond our solar system. The Very Energetic Radiation Imaging Telescope Array System (VERITAS) Collaboration joined this program in 2018 and provides the capability to search for one particular technosignature: optical pulses of a few nanoseconds in duration detectable over interstellar distances. We report here on the analysis and results of dedicated VERITAS observations of Breakthrough Listen targets conducted in 2019 and 2020 and of archival VERITAS data collected since 2012. Thirty hours of dedicated observations of 136 targets and 249 archival observations of 140 targets were analyzed and did not reveal any signals consistent with a technosignature. The results are used to place limits on the fraction of stars hosting transmitting civilizations. We also discuss the minimum pulse sensitivity of our observations and present VERITAS observations of CALIOP: a space-based pulsed laser on board the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations. The detection of these pulses with VERITAS, using the analysis techniques developed for our technosignature search, allows a test of our analysis efficiency and serves as an important proof of principle. 

Abstract While the sources of the diffuse astrophysical neutrino flux detected by the IceCube Neutrino Observatory are still largely unknown, one of the promising methods to improve our understanding of them is investigating the potential temporal and spatial correlations between neutrino alerts and the electromagnetic radiation from blazars. We report on the multiwavelength target-of-opportunity observations of the blazar B3 2247+381, taken in response to an IceCube multiplet alert for a cluster of muon neutrino events compatible with the source location between 2022 May 20 and 2022 November 10. B3 2247+381 was not detected with VERITAS during this time period. The source was found to be in a low-flux state in the optical, ultraviolet, and gamma-ray bands for the time interval corresponding to the neutrino event, but was detected in the hard X-ray band with NuSTAR during this period. We find the multiwavelength spectral energy distribution is described well using a simple one-zone leptonic synchrotron self-Compton radiation model. Moreover, assuming the neutrinos originate from hadronic processes within the jet, the neutrino flux would be accompanied by a photon flux from the cascade emission, and the integrated photon flux required in such a case would significantly exceed the total multiwavelength fluxes and the VERITAS upper limits presented here. The lack of flaring activity observed with VERITAS, combined with the low multiwavelength flux levels, as well as the significance of the neutrino excess being at a 3σlevel (uncorrected for trials), makes B3 2247+381 an unlikely source of the IceCube multiplet. We conclude that the neutrino excess is likely a background fluctuation. 

Abstract The ground-based gamma-ray observatory Very Energetic Radiation Imaging Telescope Array System (VERITAS, https://veritas.sao.arizona.edu/ ) is sensitive to photons of astrophysical origin with energies in the range between ≈85 GeV and ≈30 TeV. The instrument consists of four 12 m diameter imaging Cherenkov telescopes operating at the Fred Lawrence Whipple Observatory in southern Arizona. VERITAS started four-telescope operations in 2007 and collects about 1100 hr of good-weather data per year. The VERITAS collaboration has published over 100 journal articles since 2008 reporting on gamma-ray observations of a large variety of objects: Galactic sources like supernova remnants, pulsar wind nebulae, and binary systems; extragalactic sources like star-forming galaxies, dwarf-spheroidal galaxies, and highly variable active galactic nuclei. This note presents VTSCat: the catalog of high-level data products from all VERITAS publications. 

Abstract HESS J0632+057 belongs to a rare subclass of binary systems that emit gamma rays above 100 GeV. It stands out for its distinctive high-energy light curve, which features a sharp “primary” peak and broader “secondary” peak. We present the results of contemporaneous observations by NuSTAR and VERITAS during the secondary peak between 2019 December and 2020 February, when the orbital phase ( ϕ ) is between 0.55 and 0.75. NuSTAR detected X-ray spectral evolution, while VERITAS detected TeV emission. We fit a leptonic wind-collision model to the multiwavelength spectra data obtained over the four NuSTAR and VERITAS observations, constraining the pulsar spin-down luminosity and the magnetization parameter at the shock. Despite long-term monitoring of the source from 2019 October to 2020 March, the MDM observatory did not detect significant variation in H α and H β line equivalent widths, an expected signature of Be-disk interaction with the pulsar. Furthermore, fitting folded Swift-XRT light-curve data with an intrabinary shock model constrained the orbital parameters, suggesting two orbital phases (at ϕ D = 0.13 and 0.37), where the pulsar crosses the Be-disk, as well as phases for the periastron ( ϕ 0 = 0.30) and inferior conjunction ( ϕ IFC = 0.75). The broadband X-ray spectra with Swift-XRT and NuSTAR allowed us to measure a higher neutral hydrogen column density at one of the predicted disk-passing phases.